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Abstract 


A computer code has been developed that uses an implicit finite- 
difference technique to solve nonsimilar , axisymmetric boundary 
layer equations for both laminar and turbulent flow. The code can 
treat ideal gases f air in chemical equilibrium , and carbon tetraflu- 
oride (CF4), which is a useful test gas for hypersonic blunt-body 
simulations. This is the only known boundary layer code that can 
treat CF4. Comparisons with experimental data have demonstrated 
that accurate solutions are obtained. The method should prove use- 
ful as an analysis tool for comparing calculations with wind tunnel 
experiments and for making calculations about flight vehicles where 
equilibrium air chemistry assumptions are valid. 


Introduction 

The implicit finite-difference technique that is 
well established for nonsimilar, axisymmetric bound- 
ary layer equations involving laminar flow (see refs. 1 
and 2) has been extended to include turbulent flow 
of an ideal gas (ref. 3) and reacting gas mixtures 
in chemical equilibrium (ref. 4). Blottner (ref. 5) 
has also shown that the same technique can be ap- 
plied to flows where the fluid is in chemical nonequi- 
librium. In addition, Mayne and Adams (ref. 6) 
and Anderson, Moss, and Sutton (ref. 7) have used 
this solution technique to compute flows that include 
streamline swallowing. Thus, the implicit finite- 
difference technique provides accurate and efficient 
boundary layer solutions and is useful for analyz- 
ing boundary layer flows on many different types of 
vehicles. 

A computer code is described here that uses an 
implicit finite-difference technique to solve nonsimi- 
lar, axisymmetric boundary layer equations for both 
laminar and turbulent flow. This code has been used 
extensively to make rapid engineering and design cal- 
culations, but it can also be used to compute the 
flow over a variety of three-dimensional entry vehi- 
cles through the axisymmetric analog (ref. 8). It can 
account for ideal gases, air in chemical equilibrium, 
and carbon tetrafluoride (CF4), which is a useful test 
gas for hypersonic blunt-body simulations. This is 
the only known boundary layer code that can deal 
with CF4, which is important because of the large 
number of CF4 tests that require this unique capa- 
bility to study the effects of variable specific heats (7) 
on aerothermodynamic results. 

The purpose of this paper is to document the 
code, give a detailed description of the method of 
solution, and present comparisons of calculations 
with data from experiments. These comparisons are 
presented over a wide range of body shapes, test 


media, and local flow conditions to demonstrate the 
accuracy of the approach. 

Symbols 


An 

parameter defined by equa- 
tion (44) 

A * 

constant in equation (63) 

B n 

parameter defined by equa- 
tion (45) 

C n 

parameter defined by equa- 
tion (46) 

Op 

specific heat at constant pressure, 
ft 2 /sec 2 -°R 

Ov 

specific heat at constant volume, 
ft 2 /sec 2 -°R 

Di 

parameter defined by equa- 
tion (40) 

d 2 

parameter defined by equa- 
tion (41) 


parameter defined by equa- 
tion (51) 

e n 

parameter defined by equa- 
tion (52) 

f 

function defined by equation (12) 

9 

total enthalpy ratio, H/H e 

H 

total enthalpy, ft 2 /sec 2 

h 

static enthalpy, ft 2 /sec 2 

3 

indicator, j = 0 for two- 
dimensional flow and j — 1 for 
axisymmetric flow 

K 

parameter defined by equa- 
tion (42) 


I 


k 

thermal conductivity, 
Btu/ft 2 -sec-°R 

k :* 

constant in equation (63) 

L 

axial length, ft 

l 

parameter defined by equa- 
tion (21) 

l* 

parameter defined by equa- 
tion (20) 

l** 

parameter defined by equa- 
tion (26) 

M 

Mach number 

m , n 

grid point across boundary layer 
(see fig. 2) 

N Le 

Lewis number 

Np r 

frozen molecular Prandtl number 

N Pr,t 

turbulent Prandtl number 

N Re 

unit Reynolds number, ft" 1 

P 

pressure, lbf/ft 2 

Q,w 

general quantities defined by 
equation (16) 

Q 

convective heating rate, 
Btu/ft 2 -sec 

R 

gas constant, ft 2 /sec 2 -°R 

Rn 

parameter defined by equa- 
tion (47) 

r N 

nose radius, ft 

r,z 

coordinates shown in figure 1, ft 

r a 

radius of axisymmetric body, ft 

S 

Sutherland constant (see 
eq. (62)), °R 

T 

static temperature, °R 

T r 

reference temperature in 
Sutherland equation (62), °R 

U e 

boundary layer edge velocity, 
ft / sec 

U, V 

boundary layer velocity compo- 
nents (see fig. 1), ft/sec 

V 

normal velocity defined by 
equation (5), ft/sec 

Voo 

free-stream velocity, ft /sec 


%,y 

boundary layer coordinates 
tangent and normal to the 
surface (see fig. 1), ft 

z 

axial distance (see fig. 1), ft 

a 

angle of attack, deg 


coefficients in equation (28) 

P 

velocity gradient parameter 
defined by equation (23) 

7 

ratio of specific heats, Cp/cy 

7 

intermittency factor defined by 
equation (66) 

6 

boundary layer thickness based 
on velocity, ft 

6* 

boundary layer displacement 
thickness defined by equa- 
tion (58), ft 

€ 

turbulent eddy viscosity, 
slug/ft-sec 

0 

parameter defined by equa- 
tion (22) 


boundary layer momentum thick- 
ness defined by equation (60), ft 

A 

constant in equation (64) 


molecular viscosity, slug/ft-sec 


transformed boundary layer coor- 
dinates defined by equations (6) 
and (7) 

p 

density, slugs/ft 3 

T 

shear stress, lb/ft 2 


stream function defined by 
equation (12) 

UJ 

parameter defined by equa- 
tion (27) 

Subscripts: 


e 

boundary layer edge 

inc 

incompressible 

s 

stagnation point 

w 

wall 

oo 

free stream 

Superscripts: 


n 

mean flow quantity 

( / 

derivative with respect to r\ or 
fluctuating values 
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Method 

This section presents a description of the tech- 
nique used to solve the steady, compressible bound- 
ary layer equations for laminar, transitional, or tur- 
bulent flow over an axisymmetric body. The physical 
coordinate system used in the analysis is shown in 
figure 1. 


Boundary Layer Equations in Physical 
Coordinates 

The partial differential equations for a compress- 
ible, turbulent boundary layer, which express the 
conservation of mass, momentum, and energy, can be 
written as follows when the molecular and turbulent 
Lewis numbers are equal to 1 (see ref. 9): 


d [purl^j d (pVrl'j 


dx 


dy 


du du dp d 

m dx + ,,V dy - ~dx + dy 

§ 2=0 

dy 


, \ 9u 


dJL V — - — J (\ c N Pr \ ^ 

dx + P dy dy N Pr ^ + [i N Prt J dy 

(N P r -l\ e (N P rt~ l\l du 1 

+ M [\ *pr V N w )\ u d y) 


( 1 ) 

( 2 ) 

(3) 


(4) 


where a new normal velocity has been defined as 


_ . _ p'v' 

V = v -\ — — 


(5) 


eq. (3)) is required for the boundary layer solution 
and, along with the edge velocity U e , must be ob- 
tained from a separate inviscid solution. (See ref. 10.) 
For the types of flows considered here, the total en- 
thalpy H e is constant throughout the inviscid flow 
region. 


Boundary Layer Equations in Transformed 
Coordinates 

For help with the numerical integration of the 
governing equations, the physical coordinate system 
(x,y) can be transformed to a transformed coordi- 
nate system (£, 77). The well-known Levy-Lees trans- 
formation (ref. 9) removes the singularity in the gov- 
erning equations at x = 0 and reduces the growth 
of the boundary layer in the transformed coordi- 
nates as the solution proceeds downstream. The 
transformation can be written as follows: 


£ (x) = 

f X 

/ p e p e U e rl j dx 
'0 

(6) 

i(x,y) = 

PeUA rv p 
n/2 J J 'op e 

(7) 


From equations (6) and (7) the relation between 
the derivatives in the physical plane ( x , y) and the 
transformed plane (£,77) can be written as follows: 


d_ 

dx 


= PePeU e rl j 


d dy d 
d£ + dx dy 


( 8 ) 


d_ _ pUerj d_ 
dy s/2£ dy 


Now a stream function ip(x , y) can be defined that 
will satisfy the continuity equation (1); that is, 


The barred quantities in equation (5) represent mean 
flow values and the primed quantities represent fluc- 
tuating values. The j in equation (1) is equal to 0 for 
two-dimensional flow and to 1 for axisymmetric flow. 
These equations are solved subject to the boundary 
conditions at the wall, where y = 0 and 

u — 0, V = 0, and H = H w (x) 
and at the outer edge, where y — > 00 and 
u — ► U e (x) and H — > H e 

The above boundary conditions reflect an as- 
sumption of no slip or mass injection at the wall and 
a prescribed wall enthalpy. The static pressure that 
is assumed constant across the boundary layer (see 


pUr ° = dy ( 10 ) 

i> Vri ° - < u > 

If a nondimensional stream function is introduced of 
the form 

(£>»?) = ( 12 ) 

then from equations (8)— ( 12) , the following is 
obtained: 

purl=pU e rif' (13) 

and 

pl'.-,, w ^^g + ^-±g v /5/' (14) 


3 


From equation (13), the following is obtained: 

u = Uef' (15) 

If a general term of the form ^ ( W jfjj - ) is transformed 
using equation (9), the following can be obtained: 


where 


A, 

(w 

pU?rl 3 

d 

Lw m 

dy 

l dy) 


dg 

V dr,) 


where W and Q are general quantities. If these 
results are applied to equation ( 2 ), the streamwise 
momentum equation, the following result is obtained: 

I) < i7 > 

From the inviscid analysis it can be shown that, 

if the entropy at the boundary layer edge is assumed 
constant, 

d Pe _ „ tt dUe na\ 

- ~PeUe (18) 

Based on equation (18), the transformed streamwise 
momentum equation (17) can be written in the form 


rf" + 


where 


= 2 £ (/' ^ - f" 

' ,= '( 1 + J ) 


dj_ 


1 = 


PP 


PePe 

Oe 

P 


6=^ 


g= 2£ L dUe 

P u e dt 


(19) 

( 20 ) 
(21) 
(22) 

(23) 


When the energy equation (4) is transformed in 
a similar manner and the definition 


9 = 


if 

He 


(24) 


is used, the following result is obtained: 


/ r* \ „ \d( i** \ 

\N Pr ) 9 + [dr, \N Pr J 


+f 


g' + J 


(25) 


l** = l l + - 


e Np r 


(26) 


u = l 


\ P N Pr.t ) 

i ( i -^)4-^)](* /r ) <27) 


The boundary conditions for the transformed 
equations at the wall, where rj = 0 , are 

/ = 0, f = 0, and g = g w (£) 

and at the outer edge, where 77 — > 00 , 

f' — ► 1 and g —> 1 

Solution Procedure 

From the approach of Blottner (ref. 11) and Davis 
(ref. 12), equations (19) and (25) can be written in 
the general form 

d 2 w dw aw n /00 , 

-QrfT + ai + °* W + “3 + <*4 = 0 (28; 

where, for the momentum equation (19), 

W = f' 


<*1 = 7Z 


02 — 


-Pf 

l* 


03 = ¥ 


04 = 


06 

l* 

—2 U’ 


(29) 

(30) 

(31) 

(32) 


and for the energy equation (25), 
W = g 


Oi\ = 




I“/Npr 

OL 2 = 0 


(33) 

(34) 
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a 3 = 


a 4 = 


l**/N Pr 

-2C/ ; 

l**/N Pr 


(35) 

(36) 


The T) derivatives in equation (28) are replaced 
by finite-difference quotients that allow variable grid 
spacing in the 77 direction. This delineation allows 
grid points to be concentrated near the surface where 
the dependent variables change most rapidly. The 
£ derivatives are replaced by two-point backward 
differences. A typical grid system is shown in figure 2. 
The solution is assumed to be known at grid point 
m,n and is unknown at point m + l,n. Adams 
(ref. 13) defines a nomenclature that can be applied 
in the following finite-difference approximations: 


must be solved for the dependent variables f f and g 
under the restrictions discussed above. 

The system of algebraic equations represented 
by equation (43) is tridiagonal and may be solved 
efficiently by a method that Richtmyer and Morton 
(ref. 14) described. For a tridiagonal system, the 
following simple relation holds: 


Wm+\,n — EnW m +\,n+\ + e n 


(2 < n < TV — 1) (48) 


where 



R-2 - <Wm+i,i 

e 2 = — 

B 2 


(49) 

(50) 


/ 0 2 wA 

2{W n+1 + KW n -\ -(1 + K) W n ] m+1 

(37) 

E — A n 

(3 <n< N - 

1 ) 

(51) 

\ dr > 2 ) m4*l, n 

D 2 

B n + C n E n - 1 

/ dW\ 

$ 

+ 

1 

1 

1 

1 

£ 

+ 

(38) 

Rn ~ C n e n —i 

B n + C n E n -\ 

(3 < n < N - 

1 ) 

(52) 

V ® T ) ) 

Di 





/3WA _ Wm+l,n - W m ,n 

V ^ )m+l,n ^ 

where 


(39) 


D\ = {rin + 1 - In) + K 2 (v n ~ Vn-l) (40) 

D 2 = (Vn+ 1 ~ Tin) 2 + K(rj n - T)n- 1) 2 (41) 

K = Vn+1 ~ 7)71 (42) 

Vn ~ Vn-1 


After substitution of equations (37)-(39), the 
finite-difference form of equation (28) becomes 


For our problem, the boundary conditions at the 
wall (n = 1) and at the outer edge of the boundary 
layer (n = N) have been specified; thus, the values 
of W m +i f i and W m+ are known. The parameters 

E n and e n are evaluated by starting at n = 2 and 
moving out across the boundary layer to n = N — 1. 
The solution for W m + \ iTl is then obtained by starting 
at n = N — 1 and traversing back across the boundary 
layer to n = 2. 

Now, because the distribution of f f across the 
boundary layer is known, the transformed stream 
function / can be determined by numerical integra- 
tion of the equation 


A n Wm+\jn-\-\ + B n W jn +\, n -f- C n W m . i-Un—l — Rn (43) 


where 


B n — 


A = A I 
” D 2 + D X 

-2(1 + K) 


D 2 


C n = — - 
R n = — 0:3 4- 


D 1 

2 K _ K^ai 
D 2 D x 

n^V^ n 

A£ 


(44) 


+ a 2 + || (45) 


(46) 

(47) 


For the above system to be linear, the coefficients 
A n , B n , C n , and Rn must be treated as known quan- 
tities at point m, n. Thus, equation (43) represents 
a set of simultaneous linear algebraic equations that 


/(£,*?) = f f' U,V) dr) (53) 

JO 

that uses the trapezoidal rule and the notation from 
the boundary conditions that /' = 0 at 77 = 0. The 
physical grid position across the boundary layer is 
obtained from the equation 


V = 


PeU e ri 


^ r'Pe 

{ JO P 


dg 


(54) 


which is also integrated using the trapezoidal rule. 

In this paper, the momentum and energy equa- 
tions have been decoupled and linearized; thus, the 
solution at each station must be iterated to re- 
move the restrictions. To linearize the equations, all 
quantities on the right side of equations (29)-(32) 
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and (33)-(36) are treated as known when the a coef- 
ficients are evaluated. In the first iteration at station 
m+ 1 , these quantities are evaluated at the previous 
station m. In successive iterations, these coefficients 
are recomputed with quantities from the previous it- 
eration. This process is repeated until the difference 
in the assumed and calculated quantities is less than 
0.1 percent. Generally two to four iterations are re- 
quired for the typical step sizes used in the present 
work. Numerical experiments have shown that this 
approach yields good results. 


Boundary Layer Parameters 


After the boundary layer solution has been ob- 
tained at a given body station from the procedure 
outlined in the previous section, several boundary 
layer parameters must be determined. One quantity 
of interest is the local convective heat flux at the wall 
(y = 0 ) that, for a flow with molecular and turbu- 
lent Lewis numbers equal to 1 , is given by (see the 
appendix) 


( /i dh\ 

Qw ~ \Np r dy) w 

where Np r is the frozen Prandtl number. 


(55) 


In the transformed coordinate system f , 77 , equa- 
tion (55) can be written as 


Qw — 


IwPeHeUeHeri 

N Pr,wV 2 ? 


Quo 


(56) 


Similarly, the local shear stress at the wall is given 
by the equation 


(du\ IwPePeU^ri c „ 

r - = ' to Ui =— L (57) 

The derivatives appearing in equations (56) 
and (57) ( g f and f") are evaluated from a three-point 
forward difference at the wall. The displacement 
thickness, which is defined in physical coordinates 
x, y by the equation 



can be written in transformed coordinates £, 77 as 


is given in f , 77 as 


0 


m 


PeUert, 


f (1 - /') d V 


(61) 


The integrals appearing in equations (59) and (61) 
are evaluated from the trapezoidal rule. 


Transport Properties 

For an ideal gas and CF4, the molecular viscos- 
ity /i is calculated using the well-known Sutherland 
formula 



where appropriate values are substituted for the ref- 
erence temperature T r , the reference viscosity /x r , 
and the Sutherland constant S. For air in chemi- 
cal equilibrium, the molecular viscosity is obtained 
from the data of Peng and Pindroh (ref. 15) and a 
table look-up procedure. 

The turbulent eddy viscosity e is calculated with 
a two-layer model. The inner region is based on a 
modification of Van Driest’s analysis (ref. 16) intro- 
duced by Patankar and Spalding (ref. 17), and the 
outer region is based on the Clauser model (ref. 18) 
as applied by Harris (ref. 3). 

The turbulent eddy viscosity for the inner region i 
is given by the equation 



where k* and A * are constants assumed to be 
approximately 0.4 and 26 in the present analysis. 

The turbulent eddy viscosity for the outer region o 
is given by the equation 



where A is a constant assumed to be 0.0168, <5* nc is 
the incompressible displacement thickness, 


6 * 


Pe^e^o ® 


(0 - /') dp 


(59) 



(65) 


Similarly, the momentum thickness, which is defined 
in x, y as 


and 7 is KlebanofTs transverse intermittency factor 
(ref. 18) defined as 


6 


rn 



pu 

PeU e 



(60) 


_ 1 - erf {5 [(y/6) — 0.78]} 

7 = 2 — 


( 66 ) 
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which can be calculated approximately by the equa- 
tion (ref. 4) 


1 

1 + 5.5 (y/Sf 


(67) 


The boundary separating the inner and outer 

regions is defined as the location where 


GHil 

For ideal gases and CF4, the molecular Prandtl 
number is assumed constant; however, for air in 
chemical equilibrium, the molecular Prandtl num- 
ber is obtained from the data of Peng and Pindroh 
(ref. 15) and a table look-up procedure. The tur- 
bulent Prandtl number is assumed constant; for all 
results presented in this paper, the turbulent Prandtl 
number is taken as 0.9. 


Thermodynamic Properties 

For an ideal gas, the thermodynamic properties 
are obtained from the ideal gas equation of state 

p = „(^)j. (69) 

based on the appropriate values for the gas con- 
stant R and the specific heat Cp. For CF4 the thermo- 
dynamic properties are obtained from the curve fits 
of Sutton (ref. 19); however, for air in chemical equi- 
librium, the thermodynamic properties are obtained 
from Hansen (ref. 20) and a table look-up procedure. 

Results and Discussion 

In this section, results from the present method 
are compared with experimental data. Two gases 
are considered: air (either ideal gas or equilibrium 
chemistry options) and CF4. Comparisons are shown 
for laminar and turbulent flow. For all results, the 
boundary layer edge conditions were obtained from 
the in viscid code of reference 10, except where noted 
otherwise. 


Ideal Gas 

Surface heat transfer distributions for flow over 
a sphere are presented in figure 3 along with ex- 
perimental data measured in the Langley 31-Inch 
Mach 10 Tunnel (provided by John Micol of the 
Langley Space Systems Division). The results shown 
in figure 3(a) are for a Reynolds number of 2.4 x 10 5 
per foot and those in figure 3(b) are for 5.4 x 10^ per 
foot. In both cases the flow is laminar and the air 
was assumed to behave as an ideal gas with 7 = 1.4. 


For the computed results shown in figure 3, a grid 
of 101 points across the boundary layer was used 
with the values of rj max and K set to 7.0 and 1.04, 
respectively. This method produced a value of A77 
at the wall equal to 0.005656. For both cases the 
computed results agree well with the experimental 
data (i.e., within ±10 percent). Similar computa- 
tions were performed that used 201 points across the 
boundary layer (results not shown) to check the ad- 
equacy of the grid, but that approach was found to 
have very little effect on the computed heating rates. 

The next case considered is the flow over a spher- 
ically blunted 8° half-angle cone at M = 5 and 
a = 0°. The results for this case are presented in fig- 
ures 4(a) and 4(b) for Nft e = 2.09 x 10 c and 19.8 x 10 6 
per foot, respectively. The experimental data for this 
case are from reference 21. For the lower Reynolds 
number (fig. 4(a)), the flow is completely laminar, 
and, except for just downstream of the stagnation 
point, the present calculations are in excellent agree- 
ment with the data. At the higher Reynolds number 
(fig. 4(b)), transition starts on the nose of the body 
near x/r yy ~ 0.3, and the flow is fully turbulent near 
x/r/v ~ 0.8. This case was modeled in the present 
calculations based on a finite transition region start- 
ing at x/r yv = 0.3 that followed the approach of ref- 
erence 22. The results generally agree well with the 
experimental data. 

The next comparison, presented in figure 5, is for 
flow over a spherically blunted, 70° half-angle cone 
at Moo = 9.86, Nr 6 = 0.55 x 10 6 per foot, and 
a = 0°. The experimental data were obtained in the 
Langley 31-Inch Mach 10 Tunnel by Charles Miller 
and Ray Midden, also of the Langley Space Systems 
Division. The data are the same as those used in 
reference 23. This configuration had a sharp cor- 
ner at the trailing edge. Because of the free-stream 
conditions and the large half-angle cone, the entire 
shock layer is subsonic, with the sonic line attaching 
near the corner in a strong expansion region. The 
flow field was computed over this body based on the 
time-dependent method of reference 10 and a geom- 
etry comprising a spherically blunted 70° half-angle 
cone forebody with a trailing 40° cone frustum that 
caused the sonic line to attach at the corner. As 
the flow moves away from the stagnation point, it 
expands and the heating decreases as would be ex- 
pected. However, as the flow approaches the trailing 
edge, the heating increases rapidly, reaching a peak 
only slightly below the stagnation value. This rapid 
rise in heating near the trailing edge is the result 
of the strong expansion in this region (i.e., large in- 
crease in the velocity gradient). Some scatter is evi- 
dent in the experimental data around the stagnation 
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point, but the calculations generally yield very good 
agreement with the data over the remainder of the 
body. 

The next case considered is the turbulent flow of 
ideal gas over a flat plate. For this case the calcu- 
lated turbulent velocity profiles are compared with 
experimental data from reference 24 in figures 6(a)- 
6(c) for Mach numbers of 1.982, 3.701, and 4.554, 
respectively. The air was assumed to be an ideal gas 
with 7 = 1.4. For the computed results shown in fig- 
ure 6, a grid of 201 points across the boundary layer 
was used with the values of r] max and K set to 100 
and 1.09, respectively. This setting produced a value 
of A77 at the wall equal to 0.2944 x 10 -6 . The in- 
crease in the number of grid points was required to 
resolve the large gradients in the turbulent flow veloc- 
ity profile near the surface. For each Mach number, 
the calculated results are in excellent agreement with 
the experimental data (i.e., within ±5 percent). 

Through numerical experimentation, a grid of 
101 points with rj max = 7.0 and K = 1.04 has been 
found to provide good heat transfer results for most 
laminar flow cases; for most turbulent flows, a grid 
of 101 to 201 points with r] max = 100 and K = 1.09 
appears to provide good heat transfer results. Thus, 
these grids appear to be adequate for most problems. 

The computed results obtained thus far for per- 
fect gases and laminar and turbulent flows agree well 
with the data from experiments and indicate that the 
boundary layer code is useful for these conditions. 
In the next two sections, applications for other gas 
models are considered. 

Air in Chemical Equilibrium 

Calculated heat transfer rate distributions for 
laminar flow over a 45° spherically blunted cone with 
rj\[ = 1.0 ft are presented in figure 7 and compared 
with viscous shock layer (VSL) calculations obtained 
by Gupta, Lee, and Zoby using the method described 
in reference 25. The free-stream conditions corre- 
spond to a flight Mach number of 15 at an altitude 
of 80 000 ft. Three solutions with the present method 
are presented based on three different assumptions to 
obtain the edge conditions. The first case, labeled 
blunt cone , follows a classical boundary layer ap- 
proach and uses properties along the surface stream- 
line obtained from an inviscid, blunt-cone flow-field 
solution. Because the flow along an inviscid surface 
streamline has passed through a normal shock wave, 
the entropy along this streamline is constant and has 
a high value equal to that downstream of a normal 
shock. For this case the present results are slightly 


lower than the VSL results, with the difference in- 
creasing downstream. 

The second case, labeled sharp cone , also follows 
a classical approach but uses an inviscid, sharp-cone 
flow-field solution. All inviscid flow for this case 
has passed through a relatively weak, constant-angle 
shock wave; thus, the entropy along the inviscid 
surface streamline is constant and has a relatively 
low value (much lower than for the blunt-cone case). 
Results for this case, shown only for x/r jy > 2.0, 
are higher than the VSL and blunt-cone results. 
This approach represents an upper limit for the 
laminar heating that will be reached far downstream 
of the stagnation point when all variable-entropy 
inviscid flow has been entrained in the boundary 
layer. This limit has obviously not been reached for 
these conditions. 

The third case, labeled variable entropy , was com- 
puted by interpolation of the edge properties from 
the inviscid flow field at a distance equal to the 
boundary layer thickness away from the wall. The 
entropy at the edge of the boundary layer for this case 
is variable, starting at normal shock value at the stag- 
nation point and decreasing downstream. The heat- 
ing for this case is almost identical to the heating for 
the blunt-cone case for 0 > x/rjy < 2.0 then departs 
and follows approximately the same level as did the 
VSL results. This progression clearly demonstrates 
the influence of variable-entropy edge conditions on 
the downstream heating. The stagnation-point heat- 
ing is slightly lower (approximately 10 percent) than 
for the VSL. This difference may be attributed to 
differences in the thermodynamic or transport prop- 
erties in the two codes. 

Calculated surface heating rates from the present 
code are compared with experimental Reentry F 
flight data (ref. 26) in figure 8. The Reentry F 
configuration is a spherically blunted, 5° half- angle 
cone, 13 ft long, with an initial nose radius of 0.1 in. 
Inviscid edge conditions for these calculations were 
obtained from equilibrium-air, sharp-cone solutions 
(ref. 27). The comparisons presented in figure 8(a) 
are for 120000 ft and a Mach number of 19.25. 
The flow for this case is laminar and the calculated 
heating agrees reasonably well with the experimental 
data. 

Similar comparisons are presented in figure 8(b) 
for 80000 ft, when the flow starts to become turbu- 
lent at a 2/L location near 0.625. The calculations 
started with laminar flow at the nose and initiated 
“instantaneous” transition at z/L = 0.625. Both 
the laminar and the turbulent calculations agree well 
with the experimental data. 
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CF 4 Gas 

Surface heat transfer distributions for flow over 
a sphere are presented in figure 9 along with 
experimental data measured in the Langley 20-Inch 
Mach 6 Tunnel (provided by Charles Miller of the 
Langley Space Systems Division). The CF4 gas has 
a relatively low “effective 7” and is used to test blunt 
reentry configurations because it simulates the high- 
density ratio across the strong shock waves on these 
configurations in hypersonic flight. The flow is lami- 
nar and the predicted values agree well with the ex- 
perimental data. This case illustrates the good ap- 
plication of the CF4 chemistry package (from ref. 19) 
in the present code. 

Concluding Remarks 

An implicit finite-difference method has been used 
to obtain steady-flow solutions for axisymmetric, 
laminar, and turbulent boundary layer flow over sev- 
eral bodies and flow conditions. Results have been 


presented for laminar and turbulent flow of ideal 
gas, laminar flow of air in chemical equilibrium, and 
laminar flow in carbon tetrafluoride (CF4). Most of 
the results were based on conditions of constant en- 
tropy at the boundary layer edge obtained from an 
inviscid flow-field solution at the surface. 

Comparisons with experimental data have demon- 
strated that accurate solutions can be obtained with 
this approach. The method should prove useful as an 
analysis tool for comparisons with experimental wind 
tunnel data and for calculations of flight conditions 
wherein the assumptions of equilibrium air chemistry 
and constant entropy at the boundary layer edge are 
valid. In addition, the method can be extended to 
include variable boundary layer edge entropy effects. 


NASA Langley Research Center 
Hampton, VA 23681-0001 
October 21, 1992 
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Appendix 

Surface Heat Transfer Rate Witli 

Enthalpy Gradient 

In general, for the laminar flow of a reacting gas 
with binary diffusion, the heat transfer to the surface 
can be written in the form 

*■“(*!£ + i >p T, h ‘i%) w < AI > 

where the first term on the right represents conduc- 
tion and the second term represents diffusion. The 
symbol V is the binary diffusion coefficient, h ? ; is the 
enthalpy of the zth species, and c 2 is the mass fraction 
of the ith species. 

Now, for a reacting gas, the enthalpy is given by 


Thus, 

dT _ 1 dh 1 dc i 

dy Cp dy Cp ^ 1 dy 

Substituting this result into equation (Al) ob- 
tains the following equation for the heat transfer: 



When the Lewis number Np e , defined as 


N U = ^ (AS) 

is introduced, the surface heat transfer rate can be 
expressed as 



h = T, c * h * 

Thus, dh /dy can be written as 


dh ^ dh i dc i 


Now, 


dh, 


dT 


dy Cp ' dy 


(A2) 


(A3) 


(A4) 


When this result is substituted into equation (A3), 
the following is obtained: 


ah 

dy 


y c i°p, 


dy 2-^ ' dy 


dT sr ' , dc i 
dy ' Z^ hi dy 


dy 


(AS) 



For Np e = 1, the second term in equation (A9) is 
zero and the surface heat transfer rate becomes 


(k dh\ 

qw ~ Up dy) w 


(A10) 


or 



(All) 


where Np r is the Prandtl number. Similar results 
can be obtained for turbulent flow. 
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Figure 2. Typical grid system. 
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(a) N Re = 2.4 x 10 5 per foot; = 9.55; = 0.6235 lbf/ft 2 ; = 97.2°R. 



(b) N Re = 5.4 x 10 5 per foot; Moo = 9.74; Poo = 1.362 lbf/ft 2 ; = 95.7°R. 


Figure 3. Heat transfer rates on sphere in air. 7 = 1.4; r/y = 0.16667 ft; T w = 540°R. 



(a) N Re = 2.09 x 10 6 per foot, p^ = 34.56 lbf/ft 2 ; = 207.5°R. 


i w’ 


Btu/ft 



(b) N Re = 19.8 x 10 6 per foot. Poo = 318.43 lbf/ft 2 ; = 200°R. 

Figure 4. Heating on 8° sphere cone at a = 0°. = 5; 7 = 1.4; r/y = 0.2083 ft. 


q 


w' 


Btu/ft 



Figure 5. Heating-rate distribution on 70° half-angle sphere cone at a = 0°. M 0 0 = 9.86; poo = 1.28 lbf/ft 2 ; 
Too = 84°R; r N = 0.08333 ft; T w = 550°R. 
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(a) Moo = 1-982. poo = 259.97 lbf/ft 2 ; Too = 305.2°R; T w = 452.0°R. 



(b) Moo - 3.701. poo = 28.362 lbf/ft 2 ; = 150.0°R; T w = 401.2°R. 


Figure 6. Boundary layer on flat-plate velocity profile at x = 1.79 ft. 



(c) Moo = 4.554. Poo = 26.268 lbf/ft 2 ; = 107.6°R; T w = 374.7°R. 

Figure 6. Concluded. 



Figure 7. Heating-rate distribution on 45° half-angle sphere cone at a = 0°. Moo = 15; Poo = 58.5 lbf/ft 2 ; 
T w = 398°R; r N = 1.0 ft; T w = 1800°R. 
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(a) Altitude = 120000 ft. = 19.25; a = 0°. 



(b) Altitude = 80000 ft. = 19.97; a = -0.15°. 
Figure 8. Heating-rate distribution along cone for Reentry F. 


q , Btu/ft 2 -sec 



Figure 9. Heat transfer rate on sphere in CF 4 . Moo 
T w = 579. 6 ° R; r N = 0.1667 ft. 


6.237; poo = 3.9014 lbf/ft 2 ; 


303. 6°R 
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